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INTRODUCTION 


The  precipitation  of  magnetospheric  electrons  by  man-made  waves 
has  been  discussed  in  the  literature  for  a  number  of  years.  Helliwell 
et  al.  (1975)  used  whistler  observations  and  band  structure  in  chorus  to 
demonstrate  preferred  magnetospheric  wave  frequencies  which  were 
relaced  to  power  line  harmonics.  Bullough  et  al.  (1976)  used  satellite 
observations  of  VLF  wave  intensities  from  ground  stations  to  infer 
amplification  of  the  waves,  implicating  electrons  as  the  energy  source. 

Vampola  (1977)  presented  observations  of  precipitating  electrons  in  the 
slot  region  of  the  magnetosphere  which  were  grossly  correlated  in  longi¬ 
tude  with  the  location  of  powerful  ground-based  VLF  transmitters. 

Vampola  and  Kuck  (1978)  used  observations  of  inner  zone  electrons  in 
the  drift  loss  cone  with  a  trace-back  technique  to  locate  the  actual  trans¬ 
mitters  which  were  causing  the  precipitation.  That  particular  paper  was 
the  first  actual  correlation  of  individual  precipitation  events  with  a  single 
transmitter.  Koons  et  al.(1981)  used  structures  in  the  inner  zone  drift 
loss  cone  electron  distributions  and  synoptic  VLF  data  to  correlate  elec¬ 
tron  precipitation  events  with  specific  transmissions  from  the  VLF  trans¬ 
mitters  UMS  and  NWC.  Imhof  et  al.(1983)  utilized  satellite  observations 
of  precipitating  electrons  in  the  slot,  obtained  in  conjunction  with  special 
transmitting  patterns  from  a  ground-based  VLF  transmitter,  to  obtain  a 
one-to-one  correlation  between  VLF  transmissions  and  electrons  in  the 
local  bounce  loss  cone. 

While  it  is  now  generally  acknowledged  that  man-made  waves  are  instru¬ 
mental  in  precipitating  electrons  and  the  mode  of  interaction  is  accepted  to 
be  a  resonance  between  whistler -mode  waves  and  the  electron  gyrofrequency 
(both  appropriately  doppler- shifted  to  the  same  frame),  the  precise  details 
are  still  under  investigation.  Koons  et  al.  (1981)  utilized  wave-tracing  of  non- 
ducted  wave  propagation  to  an  interaction  region  near  the  equator.  Vampola 
(1977)  hypothesized  wave-particle  interactions  much  lower  on  the  field  line. 
Helliwell  et  al.  (1975)  utilized  interactions  in  the  equatorial  region,  but  with 
ducted  waves.  Koons  et  al.  (1981)  calculated  that  a  field  intensity  of  3  pT 
was  sufficient  to  produce  the  scattering  they  observed,  while  a  calculation  of 
electron  interaction  with  coherent  waves  reported  by  Inan  et  al.  (1Q82)  indicated 
that  detectable  electron  scattering  should  be  accomplished  with  as  little  as  1  pT. 


They  assumed  ducting.  Inan  et  al.  (19  <  <)  rerortea  in- situ  measurements 
near  the  equator  of  0.  1  to  0.  3  pT  levels  from  the  Siple  transmitter  before 
magneto soheric  amplification.  They  also  reported  observation  of  an  ampli¬ 
fied  wave  prior  to  its  first  passage  through  the  equatorial  region,  indicating 
that  wave-particle  interactions  could  occur  ivell  away  from  the  equator. 

Because  the  ionosphere  is  a  magnetized  plasma,  its  properties  are  aniso¬ 
tropic  for  radio  waves  reflecting  off  of  it  or  propagating  through  it.  Further¬ 
more,  it  is  generally  highly  attenuating.  As  a  result,  realistic  calculations 
of  wave  propagation  have  not  been  made  in  analytic  form.  Instead,  ray-trac¬ 
ing  techniques  are  usually  used.  The  advantage  of  ray  tracing  is  that  one  can 
describe  the  propagation  medium  in  either  continuous  or  discrete  form  and 
make  relatively  good  predictions  of  the  resultant  wave  propagation  vector  and 
amplitude  without  resorting  to  ducts.  Unfortunately,  ray-tracing  techniques 
predict  passage  through  the  ionosphere  and  into  the  magnetospheric  inter¬ 
action  region  at  low  latitudes  even  from  midlatitude  sources.  Yet,  conjugate 
studies  show  wave  propagation  between  conjugate  points.  Ducting  is  usually 
given  as  the  explanation.  Direct  propagation  requires  that  an  otherwise  re¬ 
flecting  ionospheric  cavity  be  sufficiently  transparent  to  waves  for  sufficient 
energy  to  leak  through  to  the  duct. 

The  question  of  ground-based  VLF  transmitter  wave  propagation  through 
the  ionosphere  into  the  magnetosphere  can  be  addressed  using  wave -particle 
interaction  data.  Luette  et  al.  (1977)  examined  the  longitudinal  variation  in  the 
occurrence  of  chorus  activity  on  the  basis  that  chorus  was  related  to  wave- 
particle  interactions  in  which  the  waves  were  power-line  harmonics.  The 
longitudinal  dependence  of  electron  precipitation  related  to  VLF  transmitters 
was  addressed  by  Vampola  (1977)  and  Vampola  and  Gorney  (1983),  but  only 
crude  correlations  with  the  location  of  gr.ound-based  sources  were  achieved. 

In  this  study,  we  will  examine  the  precipitation  pattern  around  a  single  VLF 
transmitter  during  a  single  event  to  address  the  question  of  wave  energy  entrance 
to  the  magnetosphere. 

OBSERVATIONS  AND  ANALYSIS 


The  approach  used  in  this  investigation  was  to  use  a  trace-back  technique 
(Luhmann  and  Vampola,  1977)  in  which  electrons  locally  observed  in  the  drift 
loss  cone  (i.  e.  ,  have  mirror  altitudes  such  that  they  will  interact  with  the 


atmosphere  prior  to  or  at  the  longitude  of  the  South  Atlantic  Anomaly  and 
hence  are  only  qua  si -trapped)  are  pitch-angle  analyzed  to  determine  the 
longitude  at  which  their  present  loss-cone  angle  corresponds  to  a  100  km 
atmospheric  loss -cone  angle.  This  identifies  the  longitude  of  origin  under 
the  assumption  that  sufficient  scattering  has  occurred  to  produce  a  loss  cone 
defined  by  the  local  field  geometry.  In  certain  regions,  most  notably  between 
30°  EL  and  180°  EL,  this  technique  provides  very  precise  answers,  provid¬ 
ing  the  statistical  accuracy  of  the  pitch-angle  data  is  high  and  the  opening 
angle  of  the  instrument  is  small. 

It  appears  that  this  type  of  study  is  not  likely  to  benefit  from  a  super¬ 
position  of  a  number  of  different  events.  For  one  thing,  different  resonance 
conditions  might  occur  from  event  to  event,  resulting  in  different  energy  and 
intensity  profiles  for  each  event.  Also,  the  same  equatorial  pitch-angle  distri¬ 
bution  will  provide  varying  local  pitch-angle  distributions,  depending  on  the 
position  of  the  satellite  on  the  field  line.  Combining  data  with  differing  inten¬ 
sities  and  varying  local  pitch-angle  widths  per  unit  equatorial  pitch-angle 
width  would  probably  degrade  the  final  composite  distributions  to  the  point 
that  they  would  be  useless  for  this  study.  Since  the  type  of  investigation  we 
are  conducting  here  would  not  benefit  from  a  superposition  of  a  number  of 
different  events,  a  single  event  with  particularly  good  intensity  in  the  precipi¬ 
tated  electron  peaks  was  selected.  The  data  were  obtained  by  a  magnetic  focus 
ing  spectrometer  on  the  S3-3  satellite  on  15  January  1977.  The  spectrometer 
consisted  of  two  separate  analyzing  chambers,  one  covering  the  energy  range 
12  keV  to  160  keV  in  five  differential  energy  channels  and  another  covering  the 
range  235  keV  to  1600  keV  in  seven  differential  energy  channels.  The  lower 
energy  unit  had  a  collimator  with  a  2°x5°  half-angle  response,  the  other  had 
a  5°x5°  half-angle  collimeter.  The  direction  of  scan  (using  the  spin- stabilized 
spacecraft  as  a  scanning  platform)  was  across  the  2°  direction.  A  two-dimen¬ 
sional  deconvolution  procedure  used  on  the  data  has  been  able  to  recover  struc¬ 
ture  in  a  pitch-angle  distribution  to  the  order  of  0.  1°  for  the  narrow  collimator 
and  0.5°  for  the  other.  Particle  distributions  were  sampled  each  0.9°  for  this 
event. 

While  the  data  from  the  lower  energy  range  portion  of  the  instrument  could 
have  provided  a  reconstructed  pitch-angle  distribution  to  a  higher  degree  of 
accuracy,  this  particular  event  was  characterized  by  about  two  orders  of  mag¬ 
nitude  greater  flux  in  the  2  35  keV  channel  than  in  any  of  the  others.  The  upper 
and  lower  limits  of  response  of  that  channel  are  320  keV  and  158  keV,  respec- 


tivslv.  Virtually  all  of  the  electrons  precipitated  in  this  event  were  in  this 
energy  range.  The  data  from  the  235  keV  channel  are  shown  in  Fig.  1. 

A  table  at  the  bottom  of  the  figure  lists  the  geographic  and  geomagnetic 
coordinates  of  the  satellite  at  the  time  the  data  were  obtained.  Note  that  this 
is  slot  data  and  is  at  a  considerably  higher  L -value  than  the  events  reported 
by  Vampoia  and  Kuck  (1978),  Koons  et  aL(1981),  and  in  a  similar  work  by 
Imhof  et  aL(1981)  in  which  electrons  in  this  energy  range  were  being  pre- 
cipitated  in  the  inner  zone  by  VLF  transmitters.  The  L-value  is  similar  to 
that  of  the  observations  of  Imhof  et  al.  (1983),  but  the  energy  is  higher.  The 
vertical  dashed  lines  in  Fig.  1  delimit  the  data  used  in  the  present  study. 

The  geometry  of  the  location  of  the  satellite  relative  to  the  atmospheric 
limit  of  stable  trapping  is  shown  in  Fig.  2.  The  solid  line  in  Fig.  2  corre¬ 
sponds  to  the  equatorial  pitch-angle,  a ,  of  a  particle  which  mirrors  at  100 
km  in  the  northern  hemisphere,  plotted  as  a  function  of  east  longitude.  For 
this  plot,  "East  Longitude"  refers  to  the  longitude  of  the  northern  103  km 
intercepc,  not  the  equatorial  value.  This  was  done  because  the  combination 
of  transmitter  frequency  and  particle  energy  insures  that  the  interaction  did 
not  take  place  at  the  equator,  but  rather  relatively  low  on  the  field  line.  In 
Fig.  2,  the  "X"  designates  the  longitude  of  the  satellite  and  the  equatorial 
pitch  angle  of  a  particle  mirroring  at  the  location  of  the  satellite.  For  elec¬ 
trons  to  be  observed  by  the  satellite,  they  would  have  to  have  an  equatorial 
pitch  angle  less  than  12.2°  or  they  would  mirror  above  the  altitude  of  the  satel¬ 
lite.  Also,  if  they  had  an  equatorial  pitch  angle  below  10.  5°,  they  would  be  in 
the  local  loss  cone;  they  would  mirror  below  100  km  locally.  From  Fig.  2, 
one  can  see  that  any  electrons  observed  by  the  S3 -3  satellite  must  have  had 
their  pitch  angles  perturbed  east  of  50°  EL  during  this  drift  period.  The  lined 
area  indicates  the  limit  of  observation  for  the  satellite  in  terms  of  equatorial 
pitch  angle  and  east  longitude.  (Basically,  this  defines  the  region  in  which  the 
interaction  which  produced  the  precipitated  electrons  of  Fig.  1  occurred.) 

The  method  of  determining  the  interaction  pattern  around  a  transmitter 
(and  presumably  the  radiation  intensity  above  the  ionosphere)  will  be  to  assume 
that  the  pitch-angle  distribution  observed  at  a  given  point  can  be  represented 
as  the  sum  of  a  number  of  distributions  in  which  the  original  distribution  is 
isotropic  outside  the  local  atmospheric  loss  cone  and  zero  inside,  at  the  point 
at  which  the  interaction  took  place,  and  that  the  isotropic  flux  is  proportiai  al 
to  the  intensity  of  the  waves  which  scattered  the  particles.  Note  that  since  the 
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Figure  1.  Electrons  in  the  drift  loss  cone  precipitated  by  a  VLF  trans 
mitter  on  15  January  1977.  The  data  between  the  vertical 
dashed  lines  was  used  in  this  study. 
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BOUNCE  LOSS  CONE 


satellite  can  observe  only  those  electrons  with  equatorial  pitch  angles  below 
12.2°,  the  original  distribution  need  not  have  been  completely  isotropic.  We 
require  only  that  the  distribution  be  isotropic  between  12.2°  and  the  angle 
corresponding  to  the  local  atmospheric  loss  cone  at  the  location  of  the  wave- 
particle  interaction.  The  assumption  of  isotropy  would  not  be  valid  were  it 
not  for  the  fact  that  all  of  the  particles  which  are  observable  by  the  satellite 
at  this  point  in  space  have  been  scattered  down  from  farther  up  the  field  line 
during  the  immediately  preceding  few  minutes  (during  the  time  it  takes  for 
the  particles  to  drift  from  the  point  of  scattering  to  the  location  of  the  satellite. 
Previous  observations  of  scattering  by  VLF  waves  are  consistent  with  isotropy 
in  the  region  of  scattering  in  and  near  the  local  bounce  loss  cone. 

The  method  of  reconstruction  of  the  pattern  of  precipitation  around  a  trans¬ 
mitter  utilizes  the  fact  that  the  local  bounce  loss  cone  angle  is  a  function  of 
east  longitude  and  that  the  local  pitch  angle  can  be  related  to  pitch  angle  at 
another  location  through  its  equatorial  pitch  angle.  Figure  3  shows  a  plot 
of  the  local  bounce  loss  cone  as  a  function  of  local  pitch  angle  and  equatorial 
pitch  angle.  The  dashed  line  represents  the  location  of  the  local  bounce  loss 
cone  at  110°  EL.  The  inset  presents  the  B  and  L  coordinates  of  the  satellite 
at  that  longitude.  The  solid  curve  is  marked  at  the  local  bounce  cone  angle 
for  various  longitudes  (again,  northern  hemisphere,  100  km  altitude  intercept). 
For  the  purpose  of  illustration,  an  idealized  sensor  with  a  large  geometric 
factor  and  a  zero  opening  angle  would  see  an  isotropic  distribution  above  77° 
pitch  angle  and  zero  below  for  an  electron  population  which  last  interacted  with 
the  atmosphere  at  60°  EL  and  was  unperturbed  since.  (In  actuality,  scattering 
in  the  residual  atmosphere  above  100  km  would  cause  a  few  particles  to  be 
present  with  lower  pitch  angle,  including  the  local  bounce  loss  cone. ) 

The  convolution  program  which  was  used  in  this  investigation  does  not 
actually  produce  an  input  pitch-angle  distribution  from  the  measured  data  points. 
The  program  uses  a  knowledge  of  the  two-dimensional  response  of  an  aperture 
with  finite  area,  length,  and  opening  angle.  It  includes  the  effects  of  rotational 
movement  of  the  aperture  during  data  acquisition  and  uses  both  elevation  and 
azimuthal  parameters  of  the  aperture  look  direction  with  respect  to  the  magnetic 
field  line.  It  assumes  the  particle  distribution  is  azimuthally  symmetric  about 
the  field  line.  Civen  an  input  pitch -angle  distribution  function,  it  calculates  the 
instrumental  response  to  that  distribution.  This  calculated  response  can  be 
compared  directly  with  the  raw  response  of  the  instrument  to  a  magneto  spheric 
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Figure  3. 


LOCAL  PITCH-ANGLE 


Plot  showing  the  100  km  atmospheric  loss  cone  as  a  function  of 
equatorial  pitch  angle  and  local  pitch  angle  for  an  observer  at 
B=.4396  g,  L=2.  51  at  110°  EL  in  the  northern  hemisphere.  The 
numbers  along  the  curve  refer  to  the  EL  at  which  electrons  with 
that  specific  local  pitch  angle  would  have  mirrored  at  100  km. 


distribution.  In  use,  one  makes  an  "educated"  guess  as  to  the  original  pitch- 
angle  distribution,  generates  an  output  response  with  the  program,  compares 
this  to  the  measured  distribution,  used  the  differences  to  assist  in  modifying 
the  inout  distribution,  and  iterates.  When  a  satisfactory  match  is  obtained, 
one  can  be  confident  that  the  actual  source  distribution  has  been  identified. 

However,  in  our  study  we  wish  to  relate  the  locally  observed  pitch-angle 
distribution  to  a  longitudinal  pattern  of  precipitation.  The  actual  pitch-angle 
distribution  at  any  location  is  not  of  interest,  since  it  is  only  an  intermediate 
steo  between  the  source  function  and  the  observed  distribution.  Hence,  the 
inout  was  actually  defined  in  terms  of  a  longitudinal  precipitation  pattern.  The 
selection  of  the  longitudinal  precipitation  patterns  was  somewhat  arbitrary. 

These  will  be  discussed  later. 

Since  small  differences  near  the  edge  of  the  loss  cone  will  be  needed  to  / 
distinguish  between  several  input  pattern  candidates,  it  is  imperative  that  very 
good  statistics  are  available  for  the  data  points  obtained  from  the  spectrometer. 
The  method  of  improving  on  the  statistics  was  to  forego  obtaining  information 
about  the  latitudinal  structure  of  the  radiation  pattern  and  concentrate  on  the 
longitudinal  structure  by  averaging  several  peaks  together  to  improve  the 
statistical  reliability  of  data  points  near  the  loss  cone  angle.  This  was  done 
by  obtaining  an  envelope  fit  to  the  peaks  shown  in  Fig.  1  between  the  dashed 
lines.  Note  that  in  Fig.  1  points  plotted  below  1  count/.  0625  sec  are  actually 
zero  courts  in  the  sample  and  are  plotted  only  to  show  the  number  of  samples 
in  which  no  counts  were  observed.  The  parameters  of  the  fit  were  used  to  nor¬ 
malize  all  of  the  other  data  points  to  the  peak  at  about  19405  UT.  Data  were 
then  averaged  into  1°  bins  between  1°  and  90°  pitch  angle.  For  this  data,  the 
soin  vector  of  the  instrument  was  very  close  to  perpendicular  to  the  local  field 
line  and  data  samples  were  obtained  at  all  angles,  including  those  along  the 
field  line.  The  normal  to  the  aperture  scanned  to  within  less  than  1°  of  the 
field  line.  The  above  binning  procedure  assumed  that  there  was  symmetry  in 
the  UDeomming  and  iowngoing  pitch-angle  distributions,  which  is  a  valid  assumpt¬ 
ion  provided  that  the  distribution  is  in  equilibrium  (i.  e.  ,  no  significant  local 
scattering  is  occurring).  The  large  number  of  samples  in  the  local  loss  cone 
which  show  no  flux  irrespective  of  whether  it  is  an  upcomming  or  downgoing 
sample  (Fig.  1)  shows  the  distribution  is  not  rapidly  changing.  By  this  means, 
sixteen  quadrants  of  pitch-angle  measurements  were  averaged  together,  giving 
a  typical  improvement  in  statistical  accuracy  of  a  factor  of  four.  The  statistical 
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(lcr)  limits  for  the  data  points  plotted  in  the  comparisons  with  convolution 
calculations  in  Figs.  4  and  6  are  less  than  the  size  of  the  symbol  down  to 
about  55°  pitch  angle  and  increases  up  to  about  1.  5  times  the  symbol  size 


at  50°. 


RESULTS  AND  DISCUSSION 

Figure  4  presents  a  comparison  of  the  averaged  electron  pitch-angle 
distribution  with  four  hypothetical  longitudinal  precipitation  patterns. 

Figure  4a  shows  the  distribution  (dotted  line)  as  it  would  be  observed  by 
the  S3-3  electron  spectrometer  if  pitch-angle  scattering  were  occurring 
uniformly  in  longitude  around  the  location  of  the  transmitter  with  the  ampli¬ 
tude  of  the  scattering  being  independent  of  longitude.  The  longitude  interval 
was  15°,  centered  on  the  transmitter  longitude,  as  discussed  in  Koons  et  al. 
(1981).  In  Fig.  4a,  the  accumulation  of  flux  near  the  local  bounce  loss  cone 
produces  a  response  that  is  clearly  higher  than  the  observed  response.  The 
calculated  and  observed  responses  are  all  normalized  at  90°  for  Figs.  4  and 
6.  Figure  4b  presents  a  comparison  in  which  the  input  distribution  is  assumed 
to  be  approximately  isotropic  outside  the  local  bounce  loss  cone  and  1%  of 
that  level  inside.  The  assumption  of  the  1%  backscatter  coefficient  is  based 
on  experimental  observations  during  strong  diffusion  precipitation  events. 

In  Fig.  4’o,  the  predicted  response  near  the  edge  of  the  loss  cone  is  even 
greater  than  that  of  Fig.  4a  and  deviates  strongly  from  the  observed  distri¬ 
bution.  For  the  next  case.  Fig.  4c,  a  delta -function  distribution  over  the  lo¬ 
cation  of  the  transmitter  was  assumed.  In  this  case,  the  local  atmospheric 
bounce  loss  cone  at  the  longitude  of  the  transmitter  shapes  the  pitch-angle 
distribution  and  it  remains  unperturbed  until  observed  by  the  satellite.  The  fit 
is  much  better  in  this  case.  In  fact,  it  is  excellent  except  for  a  critical  area 
at  the  edge  of  the  local  loss  cone  where  it  is  deficient  in  electrons.  Referring 
back  to  Fig.  3,  one  notes  that  precipitation  at  more  easterly  longitudes  results 
in  flux  at  lower  pitch  angles.  Hence,  Fig.  4c  indicates  that  some  additional 
scattering  must  occur  to  the  east  of  the  transmitting  site.  The  question  is, 
how  much  additional  scattering  should  we  add  and  how  can  we  justify  it. 

An  experiment  was  performed  (Likhter  et  al.f  1973)  in  which  a  ship  was 
stationed  at  the  conjugate  point  of  a  VLF  transmitter  and  then  a  second  ship  ran 
longitudinal  and  latitudinal  traverses  to  obtain  the  ratio  of  received  power  in 
the  transmitted  signal  at  the  two  locations.  For  the  present  study,  the  traverse 
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Comparisons  between  the  averaged  pitch-angle  distribution  of 
electrons  of  Fig.  1  (symbols)  and- the  calculated  instrumental 
response  for  a  precipitation  pattern  of  the  type  shown  in  the  c 


data  were  translated  into  longitudinal  precipitation  intensity  data,  centered 
at  the  transmitter  longitude,  and  used  as  an  input  source  function  by  the  con¬ 
volution  program.  The  result,  Fig.  4d,  is  again  an  excellent  fit  and  is  a  sig¬ 
nificant  improvement  over  that  of  Fig.  4c  in  the  critical  loss  cone'area.  Note 
that  for  both  Figs.  4c  and  4d  the  predicted  flux  above  30°  local  pitch  angle  is 
essentially  isotropic,  as  is  the  actual  measured  pitch-angle  distribution. 

Below  50°,  both  the  predicted  response  and  the  measured  response  (within 
statistics)  are  isotropic,  though  at  a  much  lower  intensity.  (The  calculation 
assumes  that  any  flux  which  appears  within  the  local  atmospheric  loss  cone 
is  isotropized  with  a  1%  reflection  coefficient.  )  The  small  deviation  between 
51°  and  54°  local  pitch  angle  in  Fig.  4d  may  be  due  to  scattering  by  the  resid¬ 
ual  atmosphere  during  the  drift  time  from  the  location  of  precipitation  to  the 
location  of  observation  or  it  could  be  due  to  some  other,  unknown,  scattering 
event  or  scattering  mechanism.  We  do  not  favor  scattering  by  the  residual 
atmosphere  as  the  origin  since  such  scattering  should  be  large-angle  scattering 
which  would  essentially  isotropize  the  scattered  flux  rather  than  causing  a  build¬ 
up  in  flux  just  at  the  edge  of  the  loss  cone.  The  origin  of  this  excess  is  more 
likely  due  to  a  process  which  changes  pitch  angles  slowly. 

Since  the  Likhter  et  al.data  extended  over  only  a  factor  of  ten  in  intensity, 
the  possibility  existed  that  transmitter  energy  at  lower  levels  was  present  at 
greater  longitudinal  distance  from  the  transmitting  site.  Since  no  experimental 
measurements  were  available  to  pursue  this  question,  and  since  additional  pre¬ 
cipitation  at  greater  longitudinal  distances  could  be  accommodated  by  the  dis¬ 
crepancy  between  the  observed  and  calculated  responses  of  Fig.  4c,  theoretical 
predictions  of  ionospheric  wave  intensities  were  used.  Figure  5  presents  cal¬ 
culations  of  field  intensity  in  the  ionosphere  for  two  different  waveguide  models 
at  a  frequency  of  15.  5  kHz  (which  is  appropriate  for  our  study).  The  theoretical 
predictions  are  from  Helliwell  (1965)  using  a  procedure  devised  by  Crary  (1961). 
Crary  also  calculated  exit  point  intensities  and  got  results  similar  to  those  of 
the  later  Likhter  et  al.(1973).  Note  that  Fig.  5  differs  qualitatively  in  profile  from 
the  Likhter  et  al.  profile  (inset.  Fig.  4d).  The  Crary  calculation  predicts  a  null 
in  the  pattern  immediately  over  the  transmitting  antenna.  The  null  is  not  present 
in  the  Crary  exit  point  profile. 

Since  the  wave-particle  interaction  is  presumably  occurring  relatively  low 
on  the  field  line  over  the  transmitting  site,  the  input  profile  (assuming  the  Crary 
calculation  is  valid)  should  provide  a  better  fit  than  the  exit  profile.  The  final 
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figure.  Fig.  6,  shows  the  comparison  between  the  Crary  profile  (smoothed 
as  shown  in  the  cartoon  inset  of  Fig.  6)  and  the  S3-3  data.  The  Crary  pro¬ 
file  provides  a  slightly  better  fit  to  the  data  than  the  Likhter  et  ah profile, 
but  the  difference  is  small- -certainly  smaller  than  the  confidence  limits  we 
can  place  on  the  individual  data  points  obtained  from  the  S3- 3  instrumenta¬ 
tion.  When  we  averaged  the  four  spin  periods  of  data  in  Fig.  1,  we  introduced 
a  potential  obscuring  factor  by  also  averaging  over  latitude  (the  S3-3  satellite 
is  in  a  polar  orbit).  The  equatorial  pitch  angle  of  the  atmospheric  loss  cone 
does  not  change  significantly  over  this  latitudinal  variation,  but  the  trans¬ 
mitter  field  intensity  is  probably  azimuthally  symmetric  and  the  earth's  mag¬ 
netic  field  is  slightly  skewed  (running  slightly  northwest  to  southeast).  The 
azimuthal  symmetry  would  produce  a  slight  artificial  enhancement  in  the  aver¬ 
aged  flux  at  the  longitude  of  the  transmitter  and  the  skewing  of  the  field  would 
shift  some  of  the  averaged  flux  to  slightly  western  longitudes  (based  on  the 
fact  that  the  flux  intensities  in  Fig.  1  are  slightly  higher  at  higher  latitude). 

Both  of  these  effects  are  small,  but  because  they  were  not  eliminated  from 
the  analysis  as  possible  obscuring  factors  and  because  the  difference  between 
the  fits  to  the  Crary  and  Likhter  profiles  is  small,  we  cannot  properly  address 
the  question  of  whether  there  is  a  null  in  the  transmitter  pattern  immediately 
over  the  transmitting  site  even  though  the  Crary  pattern  gives  a  marginally 
better  fit.  If  the  non-corrected  effects  mentioned  above  are  responsible  for 
the  slight  excess  of  flux  in  the  51°  to  54°  pitch-angle  distribution  comparisons 
of  Figs  4d  and  6,  a  similar  analysis  with  an  event  with  much  better  statistics 
in  a  single  traverse  of  the  0°  to  90°  local  pitch -angle  space  might  provide  the 
answer.  The  S3 -3  data  set  is  very  large  and  identification  of  a  precipitation 
event  is  a  slow  process.  It  is  not  known  at  this  time  whether  a  precipitation 
event  meeting  more  stringent  criteria  exists  in  the  S3-3  data  set. 

A  most  important  point  to  make  is  that  precipitation  intensity  as  a  function 
of  longitude  around  a  transmitter  has  been,  shown  to  be  very  similar  (possibly 
identical)  to  the  radiation  pattern  calculated  for  leakage  through  the  ionosphere 
and  to  received  signal  strength  measured  at  the  conjugate  point.  This  implies 
that  the  wave-particle  interaction  for  these  events  is  linear  or  is  a  quasi-linear 
process  rather  than  an  initiated  instability,  and  that  either  the  waves  are  ducted 
parallel  to  the  field  line  above  the.  ionosphere  or  the  wave -particle  interaction 
occurs  very  low  on  the  field  line.  Because  of  the  relatively  low  frequency  of 
the  waves  compared  to  the  energy  of  the  particle,  it  is  likely  that  the  interaction 
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Figure  6.  Plot  similar  to  Fig.  4,  but  using  the  Crary  (1961)  ionospheric 
field  intensity  pattern. 


is  taking  place  relatively  near  to  the  mirror  points  of  the  particles  where 
the  longitudinal  velocity  of  the  electrons  is  low. 

SUMMARY 

An  individual  electron  precipitation  event  produced  by  an  interaction 
between  the  electrons  and  waves  from  a  ground-based  VLF  transmitter 
ha3  been  studied  using  data  from  the  S3-3  magnetic  electron  spectrometer. 

High  resolution  pitch-angle  distribution  data  has  been  compared  with  the 
expected  distribution  from  model  longitudinally -dependent  wave-particle 
interaction  regions.  The  data  are  fit  well  with  interaction  regions  linearly 
related  to  the  field  intensity  predicted  above  a  transmitter  on  theoretical 
grounds  and  experimentally  measured  at  a  region  conjugate  to  a  transmitter. 
The  measured  pitch-angle  distributions  do  not  resemble  the  distributions 
expected  from  uniform  precipitation  in  a  broad  region  around  a  transmitter. 

A  delta -function  precipitation  pattern  immediately  over  a  transmitter  gives 
a  fairly  good  fit,  but  is  clearly  deficient  in  electrons  which  would  have  been 
precipitated  both  to  the  east  and  to  the  west  of  the  transmitter. 

The  excellence  of  the  comparison  between  measured  pitch-angle  distri¬ 
bution  and  the  distribution  predicted  on  the  basis  of  field  patterns  around  the 
transmitter  indicates  that  the  amount  of  precipitation  is  linearly  related  to  the 
field  intensity;  therefore,  the  wave -particle  interaction  is  a  linear  or  quasi  - 
linear  process,  as  opposed  to  a  process  in  which  an  instability  is  triggered 
by  the  wave. 

Finally,  the  energy-frequency  relationship  of  the  particles  and  waves  in¬ 
volved  in  the  interaction  indicates  that  the  interaction  region  is  low  on  the  field 
line,  where  the  magnetic  field  is  large  and  the  particle  has  small  velocity  along 
the  field  line  because  it  is  near  the  mirror  point. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems* 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 


Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of-f ield-of-vlew  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  Infrared  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications;  micro- 
wave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements,  diagnos¬ 
tics  and  radiometry,  microwave/millimeter  wave  thermonic  devices;  atomic  time 
and  frequency  standards;  antennas,  RF  systems,  electromagnetic  propagation 
phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  Infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 
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